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Atoms in an optical lattice, formed by a laser stand- 
ing wave, is an ideal system for studying quantum 
nonlinear dynamics. Operating at low temperatures 
and controlling the lattice parameters, experimental- 
ists now are able to tailor practically one-dimensional 
potentials and manipulate with internal and external 
degrees of freedom of atoms. Experimental study of 
quantum chaos has been carried out with ultracold 
atoms interacting with a periodically modulated op- 
tical lattice To suppress spontaneous emission 
and provide Hamiltonian quantum dynamics atoms 
are detuned far from the optical resonance. Adiabatic 
elimination of the excited state amplitude leads to an 
effective Hamiltonian for the external motion Q corre- 
sponding to a 3/2 degree-of- freedom classical system 
which has a mixed phase space with regular islands 
embedded in a chaotic sea. De Brogile waves of ul- 
tracold atoms have been shown to demonstrate under 
appropriate conditions the effect of dynamical local- 
ization which means quantum suppression of chaotic 
diffusion 0, Q • Decoherence due to spontaneous emis- 
sion tends to suppress this quantum effect and restore 
classical-like dynamics 

A new arena of quantum nonlinear dynamics with 
atoms in optical lattices is opened if we work near the 
optical resonance and take the internal dynamics into 
account. In the Hamiltonian approximation, when 
one neglects spontaneous emission (SE), the coupling 
of internal and external atomic degrees of freedom has 
been shown to produce a number of nonlinear effects 
in rigid (i.e. without any modulation) optical lattices: 
chaotic Rabi oscillations, chaotic atomic transport, 
dynamical fractals, and Levy flights 0, 0|- In real life 
the dynamics of atoms in near-resonant laser fields is 
not deterministic because of SE. The problem of in- 
terrelation between deterministic chaos and noise is 
rather general. Natural systems are subject to noise 
which, usually, acts continuously. If noise is practi- 
cally continuous and comparatively weak we can study 
in which way it affects chaotic deterministic evolu- 
tion of the system under consideration. Spontaneous 



emission is a kind of a shot noise which is not small 
because SE recoils may change the internal state sig- 
nificantly. In this Letter we demonstrate analitically 
and in a Monte Carlo simulation that manifestation 
of dynamical instability and Hamiltonian chaos in the 
fundamental near-resonant matter-radiation interac- 
tion can be found in the behavior of atoms moving in 
a rigid optical lattice. 

Since we study manifestation of quantum nonlinear 
effects in ballistic transport of atoms, when the av- 
erage atomic momentum is very large as compared 
with the photon momentum hkf, the transnational 
motion is described classically by Hamilton equations. 
We start with the Hamilton-Schrodingcr equations of 
motion for a two-level atom in a standing light wave 
which have been derived in Refs. @, 
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where x = kfX and p = P/hkf are normalized atomic 
center-of-mass position and momentum, u, v, and 
z are synphase and quadrature components of the 
atomic electric dipole moment and the population in- 
version, respectively. The length of the Bloch vector, 
u 2 + v 2 + z 2 = 1, is conserved. The dot denotes differ- 
entiation with respect to the normalized time t = Qt. 
The values of the normalized decay rate 7 = T/fl and 
the recoil frequency u> r = hk 2 /m a fl are chosen to be 
7 = 3.3 ■ 10 -3 and oj r = 10 -5 and correspond to a ce- 
sium atom (A a = 852.1 nm and T — 3.2 ■ 10 7 s _1 ) in a 
strong field with the Rabi frequency Vl = 10 10 s — 1 . So, 
the normalized detuning between the field and atomic 
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frequencies, A = (ojf — Lo a )/VL , is a single variable 
parameter. In Eqs. §Q Tj are random time moments 
of SE events and pj are random recoil momenta with 
the values between ±1 (ID case). In terms of the nor- 
malized time r the mean frequency of SE events is 
7(2 + l)/2. At t — Tj, the atomic variables change as 
follows: p — > p + pj, u — > 0, v — ► 0, z — > — 1. 

Equations fl} with 7 = and without the terms 
containing delta-functions describe Hamiltonian co- 
herent evolution of the internal and external degrees 
of freedom of an atom that has been shown [J] to 
be chaotic (in the sense of exponential sensitivity to 
small changes in initial conditions) in certain ranges 
of values of the parameters cu r and A and initial mo- 
menta. With comparatively small values of the initial 
atomic momentum po, atoms may wander in an op- 
tical lattice with alternating trapping in the wells of 
the optical potential and flights over its hills. It is a 
kind of a random walking that may occur without any 
modulation of the lattice parameters and/or any noise 
like SE US- 

In this work we consider only fast ballistic atoms 
which never change the direction of motion. There is 
a range of large values of initial momentum p$ where 
the maximal Lyapunov exponent A of the Hamilto- 
nian equations of motion has been computed to be 
positive Q . It means that the momentum of a ballis- 
tic atom without SE may oscillate in a deterministic 
but chaotic way around a mean value (p). The cen- 
tral question of the present study is the following. In 
which way the Lyapunov instability and Hamiltonian 
chaos, that may occur between SE events, manifest it- 
self in ballistic atomic transport which is a stochastic 
process due to SE? 

To answer the question we simulate Eqs. ([T]) by a 
Monte Carlo method (for details see Q ) and compute 
atomic trajectories in the momentum space to find 
the momentum diffusion coefficient D p as a function 
of the momentum p. The results are compared with 
the maximal Lyapunov exponent A computed with the 
Hamiltonian anologue of the set {H (without SE). More 
exactly, as a measure of Hamiltonian chaos, we com- 
pute chaos probability A = (20(A) — 1), where 0(A) is 
a Heaviside function, which is equal to for A < 0, 
1/2 for A = 0, and 1 for A > 0. The values of A in 
Fig. [T] have been computed by averaging over many 
atomic trajectories with close values of p. If A > 
with all those atoms, then A = 1, and we have Hamil- 
tonian chaos with probability 1. If A = then A = 0, 
and the motion is regular with the probability 1. The 
values in the range < A < 1 mean that chaotic 
and regular trajectories are mixed in a small range 
of values of initial momenta, and chaos probability is 
proportional to the fraction of atoms with positive As. 
Fig. [T] demonstrates a correlation between the regimes 
of chaotic (regular) Hamiltonian transport and the be- 
havior of the momentum diffusion coefficient D p for 



spontaneously emitting atoms. Beginning with those 
values of the momentum p, for which the probability 
of Hamiltonian chaos becomes smaller than 1 (see the 
lower panels in Fig.[T]), one observes an abrupt transi- 
tion to a more regular regime of motion with another 
law of decay of D p (see the upper panels in Fig. [T|). 

We stress that the atomic transport in reality is 
stochastic with all the values of p due to SE, but the 
measure of its stochasticity, D p , decays rapidly in the 
same range of the momenta where the Hamiltonian 
analogue of the system demonstrates a transition from 
chaos to order. It is more important that this differ- 
ence could be measured in real experiments and would 
provide us with direct signatures of atomic Hamilto- 
nian chaos in terms of transport characteristics which 
are more easy to measure than the Rabi oscillations. 

In what follows we estimate the diffusion coeffi- 
cient D p when the corresponding Hamiltonian ballistic 
transport is chaotic and regular. In the weak Raman- 
Nath approximation, w r p 2 /2 ^> |ucosx+Az/2|, when 
the atomic kinetic energy is not strictly a constant, 
but much larger than the potential one, the momen- 
tum fluctuations between SE are small. In Ref. 0] we 
have shown that at small detunings, |A| <C 1, the evo- 
lution of the total energy H = oj r p 2 /2 — ucosx — Az/2 
(which is a constant in the absence of dissipation) is 
a quasi-random process with sudden changes in H, 
when SE occurs, and a slight linear drift in between. 
We can treat the evolution of H as the following map- 
ping: 
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(2) 

where Hj is a value of the energy just after j-th SE, 
x(Tj), u{Tj), z(jj), and p{rj) are values of the corre- 
sponding variables just before j-th SE which are de- 
termined by the evolution between SE events. The 
last term with the averaging over a time exceeding 
the period of the Rabi oscillations is a result of an en- 
ergy drift between SE events. In general, this random 
walk is asymmetric. There is a friction force F = (p) 
which can accelerate or decelerate atoms in average. 
The measure of momentum fluctuations in an atomic 
flight (duration of which exceeds largely (Tj — Tj_x)) 
is a momentum diffusion coefficient D p which can be 
written as 



((Hj-Hj.x^-iHj-H^xY 
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where the average value of the momentum in the weak 
Raman-Nath approximation is p ~ y/2H/ui r . Using 
the largest second and fifth terms in Eq. ([2|), we can 
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FIG. 1: Correlation between the average momentum diffusion coefficient in a log-log scale (in units of h 2 k 2 Q) and 
probability of Hamiltonian chaos A in their dependencies on atomic momentum p (in units of Hkf) at A = — 0.01 and 
A = —0.0005. Dashed lines with the slope p~ 2 + const and p^ 1 + const are theoretical curves l|6]l. (0, valid in the 
regimes of Hamiltonian chaos and order, respectively. Solid line is a theoretical curve (|9]) derived to fit the numerical 
dependence. An abrupt change in the decay laws for D p occurs for those values of p where the transition from order to 
chaos takes place in the Hamiltonian dynamics. 



estimate D p as follows: 
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All the other terms in Eq.Q are small since |A| <C 1, 
\z\ ~ 1, and \u\ 3> |A|. In deriving Eq. we put 
(ucosx) ~ 0, (Tj — Tj-i) — 2/7, and (p 2 ) = 1/3. 

To estimate the value of (m 2 (tj)} in Eq. (0| we use 
the results of our theory [5| where we have shown, 
that in the absence of SE the variable u can be ap- 
proximated by a constant when atoms move between 
nodes of a standing wave (it fact, it performs shallow 
oscillations) which changes suddenly its value when 
they cross any node at cos x = 0. Spontaneous emis- 
sion results in additional jumps, u — > 0, but between 
SE events one can approximately describe the dynam- 
ics using the Hamiltonian theory In the chaotic case, 
the evolution of u can be approximated as a stochastic 
mapping 



uj r p 



(5) 



where u m are values of u after m-th node crossing 
(starting with the latest SE event), <fi m are random 
phases in the range [0, 7r]. The index m increases by 
1 just after each node crossing and jumps to zero just 
after SE event. This map is obtained from the ex- 
pression (11) in Ref. [3 in the limit \u\ « 1 which 



is valid because we have u = after any act of SE , 
and the values of u never go far away from zero due to 
small magnitudes of jumps. Between the acts of SE, 
a sequence of values of u looks like a Markov chain of 
random jumps where the next state depends only on 
the previous one. In the weak Raman-Nath approxi- 
mation, the number of node crossings between SE can 
be estimated in the average to be (M) = 2u> r p/ (jtt). 
Now we can estimate the value of u(rj) ~ um and, us- 
ing Eq.Q, get the following formula for the momen- 
tum diffusion coefficient in the regime of Hamiltonian 
chaos 



,J_ 

12 



(0) 



In Fig. [T] this function ([6]) is shown by the dashed 
lines in a log-log scale. It fits well numerical data in 
the range of atomic momenta where Hamiltonian dy- 
namics is fully chaotic, i.e. at A = 1. The formula © 
is valid in a wide range of moderately small detunings, 
but it does not work in the Hamiltonian mixing and 
regular regimes. 

At very small detunings, we may estimate D p both 
in the fully chaotic (A = 1) and regular (A = 0) 
regimes. With fast atoms at |A| <^ 1, we can 
use the strong Raman-Nath approximation (neglect- 
ing the momentum fluctuations between SE at all) 
and adopt the simple linear law of motion x ~ u) r pr. 
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In Ref. [H we have derived a formula (see Eq. (A3) 
therein) for the value of u after crossing the first node. 
So, we get for u after crossing the m-th node 



Vq cos 
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where vq and z$ are constant values of v and z at x = 
irk, k = 1,2,.... Now jumps of u are not random, and 
in the range between two SE events a trajectory looks 
like a ladder with odd and even jumps of different size 
with the total number of steps M. It can be shown 
that u(rj) ~ um — MAvq \/ tt /uj r p cos(2 / uj r p — u /4) . 
Since Vo differs after different SE events, we put {vq) ~ 
1/2. At very small detunings, the diffusion coefficient, 
corresponding to Hamiltonian regular motion (A = 0), 
is 
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Thus, we have the analytic expressions for D p in the 
regimes of Hamiltonian chaos (A = 1) and Hamilto- 
nian order (A = 0). In general case, < A < 1, we 
suppose a linear law for the momentum diffusion: 



D p ~ (1 - A)D 
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This function is shown by the solid line in the upper 
right panel in Fig. [TJ 

Let us consider a small cloud of atomic gas mov- 
ing in one direction with the mean momentum (p). 
Initial position and momentum distributions are sup- 
posed to be Gaussian with the standard deviations 
a 2 = ((x — {x)) 2 ) and a 2 = {(p - {p}) 2 )- The momen- 
tum diffusion coefficient is D p = d(a 2 ) / (2dr) . The 
temperature of atomic gas and its heating speed (in 
units K/s) are 



h 2 k 2 a 2 p 
m n ku 



dT 
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where Ek is the atomic kinetic energy (in J) in the 
center-of-mass moving frame. The heating speed is 
proportional to D p which has been shown in this Let- 
ter to demonstrate different behavior in the regimes 
of regular and chaotic Hamiltonian dynamics. 

In real experiments the measurable quantity is a lin- 
ear cloud size L = 2a x /kf (in meters). At small obser- 
vation times r <C p/F and small temperatures a p <^ 



(p), we can approximate D p by a constant for all the 
atoms in a cloud which does not change significantly 
under the action of the force F. In this approximation 
we find: a 2 x = a 2 (0) + {1/2) uj 2 a 2 (0)t 2 + (2/3) D p lo 2 t 3 . 
We have computed L with Eqs. (p} and with that for- 
mula and found a good correlation between the results 
and a strict difference in extension of atomic clouds in 
the presence and absence of Hamiltonian chaos. 

In conclusion, we have found in numerical experi- 
ments the manifestation of dynamical instability and 
Hamiltonian chaos in ballistic motion of two-level 
atoms in a near-resonance standing-wave field. The 
effect of dynamical chaos in the fundamental atom- 
light interaction is masked by random events of SE. 
Nevertheless, we proved analytically and numerically 
that, under certain conditions, there exists a clear cor- 
relation between the behavior of the momentum diffu- 
sion coefficient D p and chaos probability A. To detect 
and quantify this effect in a real experiment, we pro- 
pose to measure linear extensions L of atomic clouds 
with different values of the mean atomic momentum 
(p). We predict that beginning with those values of 
(p) , for which Hamiltonian chaos probability becomes 
to be 1, the value of L for the corresponding atomic 
clouds should increase sharply. 
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